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^sj ' Abstract We propose a modified procedure for extracting the numerical value for the 

strong coupling constant cts from the r lepton hadronic decay rate into non-strange 
particles in the vector channel. We employ the concept of the quark-hadron duality 
specifically, introducing a boundary energy squared Sp > 0, the onset of the perturba- 
tive QCD continuum in Minkowski space [Tl[2][3] . To approximate the hadronic spectral 
function in the region s > Sp, we use Analytic Perturbation Theory (APT) up to the 
fifth order. A new feature of our procedure is that it enables us to extract from the 
data simultaneously the QCD scale parameter Ajj^ and the boundary energy squared 
Sp. We carefully determine the experimental errors on these parameters which come 
from the errors on the invariant mass squared distribution. For the MS scheme cou- 
Mh| pling constant, we obtain as{mT) = 0.308 ± O.OMoxp.- We show that our numerical 

analysis is more stable against higher-order corrections than the standard one. The ex- 
tracted value for the duality point Sp is found surprisingly stable against perturbation 
theory corrections sj = 1.71 ± 0.05cxp ± 0.00th GeV .Additionally, we recalculate the 

CO ' "experimental" Adler function in the infrared region using final ALEPH results. The 

"^ , uncertainty on this function is also determined. 

r^ ■ Keywords tau lepton decay • renormalization group equation • perturbation theory 

^O ' data analysis 
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1 Introduction 



The hadronic t decays serves as an ideal laboratory for testing quantum chromody- 
namics (QCD) in a relatively low energy regime. In the past, various techniques (fixed 
Jv^ , order perturbation theory, contour improved perturbation theory, effective charge ap- 

» I ■ proach, renormalons, dispersive approach) have been devised to improve the reliability 

C^ I of the predictions of the theory for the r system. In this boundary area of the en- 

ergy, perturbative ideas are still applicable due to relatively large mass of the r lepton. 
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while non-perturbative effects are expected to be small. Usually, they are under con- 
trol within Wilson Operator Product Expansion (OPE) [4]. It is known that the main 
calculational tool in perturbative QCD (pQCD) the renormalization group improved 
perturbation theory augmented with the OPE can not be used locally in the time- 
like region even at high energy. Fortunately, this problem has been resolved in earlier 
work [5] by means of the idea of the quark-hadron duality. This enabled one to employ 
the QCD perturbation theory in Minkowski region to calculate some global (inclusive) 
quantities like r lepton decay rate. Although the quark-hadron duality cannot be justi- 
fied rigorously from the first principles, in practice this idea works good enough. Using 
the duality, an accurate description of the r lepton decay data was achieved (see the 
seminal work [B] and the literature therein). However, one should always keep in mind 
that the duality between a physical quantity and its quark-gluon perturbation theory 
representation is only approximative and thus it must inevitable be violated (see the 
review [7] and the literature therein). To identify general mechanism of possible Duality 
Violations (DVs), special QCD inspired models for the hadronic spectral functions (e.g. 
the instanton-based and resonance-based models [7] as well as the models motivated 
by the large Nc limit of the theory [B.) have been studied. In these models DVs in fact 
occur. Presumably, DVs arise due to the lack of the convergence of the OPE on the 
Minkowski axis. If this is the case, then the analytical continuation of the truncated 
OPE series from the Euclidean region to the physical axis is questionable [7]. 

In recent years, the accuracy of the measurements of the observables of the r lepton 
system has been essentially improved (for the recent results of the ALEPH collaboration 
see [9l ll0l[mil2) l. This enables one to extract the parameters of the standard model 
from r data with very high precision. Of particular interest is the numerical value of the 
strong coupling constant as . Admittedly, one of the most precise determinations of the 
strong coupling constant comes from the analysis of the r data (for most recent results 
see [12]). An independent low-energy highest-precession determination of as comes 
from lattice QCD simulations combined with experimental data for hadron masses 
[13| . These two highest-precision determinations extrapolated to the Z mass yield 

Qs(M;?) = 0.1212 ±0.0011 (r decay) (1) 

Qs(M^) = 0.1170 ±0.0012 (lattice). (2) 

Note that the agreement between these two results, with the errors quoted, is not good. 
They differ from each other by about 2.6 standard deviations. Furthermore, the lattice 
determination is closer to as{Mz) values obtained from high energy experiments. Thus, 
the reliability of the estimates from the r-lepton data has been called in question |12l 
I14II15II16U17U18U19U20U21| . The small but still significant non-perturbative effects have 
been included into analysis [181121] . On the one hand, the impact of the higher order 
terms of the OPE (neglected in the standard analyzes) has been estimated [14II15||T^ 
117] . It was confirmed that their infiuence on the extracted value of as is not small 
in the separate vector and axial vector channels. To suppress these contributions in 
the finite energy sum rule the so-called pinched weights introduced [14lll5lfT6lll7j . An 
independent estimation of possible non-perturbaive corrections to the finite energy 
sum rele (direct instantons, duality violation and tachyonic gluon mass) which cannot 
be described within OPE can be found in jl^. To estimate systematic effects from 
DVs, recently the authors of [12] have analyzed the ALEPH r data for the V±A 
spectral function using two different models of DVs. These models were previously 
considered in . It was confirmed (within this models) that DVs effects in this channel 



are completely negligible. However, this problem has been reconsidered in [20) . There 
the separate vector (V) and axial-vector (A) spectral data have been analyzed. To 
describe DVs coming from the region s > 1.1 GeV physically motivated models for 
these spectral functions have been suggested. Analyzing the r data provided by the 
ALEPH collaboration, the authors of [20] have concluded that DVs are not small. An 
additional systematic error in the value of the coupling constant coming from DVs has 
been estimated on the level Sasinq.) ~ 0.003 — 0.010. 

As is well known, in the time-like region the renormalization group (RG) invariance 
cannot be used unambiguously. Usually, the QCD corrections to the r lepton decay 
rate Rt is expressed via the contour integral of the associated Adler function multiplied 
by the known weight function. This representation is valid owing to special analyticity 
structure of the corresponding exact current-current correlation function. The Adler 
function is represented via the truncated perturbation theory series and the integral is 
taken over the circle of radius m^ {rriT stands for the r-lepton mass) in the complex 
energy squared plane [5]. One possibility is to integrate term- by-term the truncated 
perturbation theory series over the contour and then perform the RG improvement. 
This approach is referred to as fixed order perturbation theory (FOPT). Alternatively, 
one can insert the RG improved truncated series for the Adler function inside the 
contour integral and then perform the integral. This approach suggested in [22112311^ 
was termed contour improved perturbation theory (CIPT). The advantage of CIPT 
is that it enables to resume some higher-order contributions to the rate. These two 
approaches lead to differing results. The values of as extracted from r decays employing 
CIPT have always been higher. A detailed comparison of these two approaches may be 
found in recent works |25II26) . A practical review of various approaches to the r decay 
rate may be found in [11| . 

The inclusive quantity like hadronic r decay rate may be accurately described 
within pure perturbative approach, provided the DVs are small. Indeed, in the V-l-A 
channel, the nonperturbative power suppressed contributions described by the OPE 
(continued analytically to the time-like region) have been estimated to be small [511121 
118) . However, the large value of the running coupling parameter at the r lepton mass 
scale leads to the large renormalization scheme dependence of perturbative predictions. 
To reduce this dependence various resummation techniques have been developed (see, 
for example, I27II28II25] '). In [27], the V-l-A r-lepton decay data was analyzed within a 
modified extraction procedure based on the effective charge approach. The numerical 
analysis has been performed in the internal renormalization scheme of the r system 
and then the result was translated into the MS scheme using renormalization scheme 
transformation. This procedure yields smaller value for the coupling constant. Similarly, 
in )28) and [29] in calculations of the t decay rate the minimal sensitivity and effective 
charge schemes were used. In this way the reliability of the estimates for the coupling 
constant has been improved. 

A serious shortcoming of the conventional perturbation theory approximations to 
the current-current correlation functions parameterized in terms of the running cou- 
pling is that they do not obey correct analytical properties of the corresponding exact 
quantities. The analytical properties are violated due to the non-physical Landau sin- 
gularities of the perturbative running coupling which appear at small space-like mo- 
menta (for the analytical structure of the perturbative coupling beyond the one-loop 
order see ,^30,31,32,33 ). Supposedly, these singularities may deteriorate the extracted 
values of the parameters [34j- This problem does not arise within dispersive or ana- 
lytic approaches to pQCD. At present, several such approaches are being intensively 



developed [3 M 3 W 37 l |38 l [39l|40l[4T]|42 p 3 p 4 p 5 p 6 p 7l[48l|49l[50l|5^ . In works [38] and 
[39| . the T lepton decay rate has been analyzed within a simple and effective disper- 
sive technique, the Analytic Perturbation Theory (APT) (for reviews see [4UII41IB51 
I46II52] ). However, the minimal analytic QCD model (the same APT) predicts, from 
the non-strange r lepton decay data, too large value for the strong coupling constant, 
as{mr) = 0.403±0.015 [39]. The advantages and shortcomings of the three approaches 
to the r decays (FOPT, CIPT and APT) were thoroughly analyzed in gj. It should 
be noted that the APT as well as its generalized versions suggested more later |47l 
148 1149 1[50] proved to be very useful from the phenomenological point of view. A re- 
markable feature of these modified expansions is the better convergence and improved 
stability property with respect to change of the renormalization scheme. Nevertheless, 
one should keep in mind that an analytic approach based only on perturbation theory 
can not be defined unambiguously, since there is not a unique recipe for removing the 
Landau singularities from the running coupling. 

A particular problem emerges from the observation that the QCD perturbation 
theory augmented with the OPE fail to describe the detailed infrared behavior of the 
Adler function associated with the r decay rate 0]. To treat this problem a more 
general framework is required. A suitable theoretical framework was suggested in [3]. 
There the hadronic non-strange vector spectral function vi (s) |J was represented by a 
simple ansatz 

vi{s) « e{sp - s)^;'/P■(s) + 6{s - Sp)i;P'^^°(s), (3) 

where v'^ (s) is the perturbation theory approximation to the spectral function and 
Sp is the onset of perturbative continuum Q an infrared boundary in Minkowski re- 
gion above which we trust pQCD. The non-perturbative component of the spectral 
function v^^'(s) was described by a resonance based model ("the lowest meson dom- 
inance approximation to large- A'^c QCD"). Using this model the authors of [3] have 
achieved correct matching in the intermediate region between the pQCD and Chiral 
Perturbation Theory predictions for the Adler function [j. To compare the Adler func- 
tion evaluated from ^ to the experiment the authors of [5] have also constructed the 
"experimental" spectral function 

^^i'^'^P" (.) = e(.p - s)vl^Pis) + e{s - Sp)^;P'5^°(s), (4) 

where ii™P(s) is the genuine experimental part of the total "experimental" spectral 
function which is measured with high precision by ALEPH [9ll53| and OPAL [54] col- 
laborations in the range < -^s < rrir = 1.777 GeV . Formula ([4]) extends the spectral 
function beyond the range accessible in the experiment. Formulas ^ and Q provide 
practical realizations of the concept of the quark-hadron duality (see the original works 
[T][2]). The ansatz ([3]) may be considered as an alternative for the truncated OPE in 
Minkowski region. The conventional formulation of the duality may be recovered from 
formulas ((3J| or Q by taking the limit Sp — >■ and introducing the OPE contributions 
IJ. Note that, the non-perturbative corrections to the spectral function described by 
model ((3]l are essentially confined in the low energy region < s < Sp. 

^ We use the normalization of the spectral function with the naive parton prediction 

"^1, naive — ^/^' 

^ The inequality < Sp < m^ is assumed. 

^ The infrared behaviour of the Adler function was also correctly described within APT 
1451 . However, to reproduce the r data, APT requires large effective quark masses. 

* Strictly speaking this is true if the perturbation theory component of the spectral function 
^pQCD^^^ is evaluated within FOPT or APT. 



In this paper we concentrate on formula (|4|. Our aim is to utilize the total in- 
formation encoded in this representation. We recall that the authors of [3] have used 
ansatz Q to extract the numerical value for the parameter Sp from the experimental 
data. For the MS scheme scale parameter (for the three active flavours) they used the 
estimate 

^MS = (372 ± 72) MeV. (5) 

The QCD component of the spectral function, v^ (s), was determined from the 
order ©(af) approximation to the Adler function. The approximation was constructed 
in terms of the exact numeric two-loop running coupling constant, normalized at the 
scale Sp. The experimental component v^^^'{s) was reconstructed from the ALEPH 
collaboration data obtained in 1999 '53^. Note that the estimate ((5)1 is close to the 
ALEPH result for the scale parameter obtained for that time 

^MS ^ (370 ± 51) MeV. 

However, these two results for Ar-r^ should not be compared. The final result of the 
collaboration for the coupling constant corresponds to the average of the two values 
obtained within the FOPT and CIPT approaches, while authors of [3] used only FOPT. 
Furthermore, in the ALEPH analysis the estimate for 0(0^) term was also included, 
while the QCD scale parameter was extracted using the exact (numeric) four-loop 
running coupling. Using the ansatz Q the authors of [3] have derived consistency 
condition from the OPE , an equation relating the parameters Sp and Aym;- From this 
equation, with the estimate ([S]), they have found that 

Sp = (1.60 ±0.17) GeV^. (6) 

Usually, it is more convenient to compare the time-like experimental data with theory 
via the Adler function, the object determined in the space-like region "^^' ' 
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^2^ _ ^2 / 2vi(s)ds 
is + Q2) 



DiQ') = Q' P^^9^^ (7) 



for this quantity reliable approximations are constructed in pQCD, in massless |56l 
I57ll58li59.60 as well as in massive cases |60II55) . The "experimental" Adler function is 
obtained by inserting ansatz Q into integral ((Tj 

D"c^p" (0^) = i3cxp(Q^ Sp) + DpQCBiQ^ Sp), (8) 

where the experimental and perturbation theory components of the total "experimen- 
tal" Adler function are defined as 

2 ^ 2 f" 2vl''P{s)ds , 2 . 2 f°° 2«PQ^°(s)ds 

(9) 
Note that the "experimental" Adler function is not wholly experimental quantity, since 
it depends also on the theoretical component Dpqqd{Q , Sp). The latter may be cal- 
culated using different theoretical approaches. For example, one may apply FOPT or 
APT. Furthermore, the result will depend on the higher order corrections to the /3- 
function and to the Adler function. In the past years, the "experimental" Adler function 



and Q > for space-like rtioinenta 



was employed for testing various theoretical approximations to the Adler function [3l 

mm- 

In view of appearance of final ALEPH data in 2005 [9lll0) it is worthwhile to recal- 
culate the "experimental" Adler function. In this paper, we will use different strategy 
for extracting numerical values of the parameters from the data. The distinguishing 
feature of our analysis is that we will determine both parameters (^T|g and Sp) self- 
consistently. Secondly, we pay particular attention to the estimation of the experimen- 
tal errors on the parameters and Adler function. Furthermore, we will use a dispersive 
approach q 

In Sect. 2 we evaluate the perturbative component of the hadronic spectral function 
up to order 0(as) within the dispersive approach. Then, we derive a transcendental 
system of equations for the parameters Ay^ and Sp. The first equation of the system fol- 
lows from the OPE for the current-current correlation function in the limit of massless 
quarks. The second equation for the parameters is a consequence of the quark-hadron 
duality implemented by means of the ansatz ^ ; perturbation theory is used to calcu- 
late the decay rate of the r-lepton into hadrons of invariant mass larger than y^. In 
Sect. 3 we solve the system of equations for the parameters numerically. To determine 
the empirical contributions in these equations, we employ the final ALEPH data on the 
non-strange vector invariant mass squared distributions which are available in [10]. To 
test the stability of the numerical results against the QCD perturbative corrections, we 
use different approximations to the Adler function from order 0{as) to order 0{ag). 
This enables us to determine the indicative theoretical errors _27_ on the extracted 
numerical values of the parameters. Our approach, which we refer to as APT , is com- 
pared with the standard CIPT. In the most of the calculations, we use the four-loop 
running coupling. In Sect. 4, we present numerical results for the "experimental" Adler 
function obtained from the final ALEPH data. The values and associated experimental 
errors of the function are tabulated in the region Q = — 1.5 GeV. Our conclusions 
are given in Sect. 5. In Appendix A we give some practical formulas obtained from the 
explicit (series) solution to the higher order RG equation. The statistical errors on the 
parameters are carefully estimated in Appendix B. In Appendix C we present some 
required results obtained within standard CIPT. 



2 Theoretical Framework 

The main quantity of interest for following analysis is the Adler function associated with 
the vector current two-point correlator. The perturbative expansion of this function in 
the limit of vanishing quark masses reads |25] 

oo n+l 2 

D{Q ) = / 'j"(m ) / kCfij^L^ where L = In — j, (10) 

n=0 k=l 

as(At ) = "°y with as{fi ) being the strong coupling constant renormalized at the 
scale fj,. Since the Adler function is a physical quantity, it satisfies a homogenous RG 
equation. This fact enables us to choose /i — Q . Then the expansion (|10|l may be 



° The difference between our framework and APT of Shirkov and Solovtsov is clarified in 
Sect. 2. 



reexpressed as an asymptotic expansion in powers of the running coupling as (Q ) 



DKGi{Q')-j:dj^^^^^] , (11) 



'as{Q 

where dn = Cn,i and the subscript "RGI" refers to the renormalization group improved 
perturbation theory. The first two coefficients in series pi[) are universal d^ = di = 1. 
The coefficients of order a^ and af in the MS scheme have been calculated about thirty 
and fifteen years ago [56II57I[58| . Recently, the authors of [59) have calculated the coef- 
ficient d4 in the case of massless quarks by using powerful computational techniques. 
The known higher order coefficients in the MS scheme for n f = 3 quark flavours take 
values d2 ~ 1.6398, d^ ~ 6.3710 and d^ ~ 49.0757. 

In practice the series (|lip should be truncated. The obtained approximations to the 
Adler function do not obey correct cut-plane analyticity properties of the exact function 
because of the non-physical "Landau singularities" which present in the perturbative 
running coupling. The exact Adler function D{z) {z — Q — —q ) is known to be 
analytic except the cut running along the negative real axis. This fact enables us to 
calculate the hadronic non-strange vector spectral function from the Adler function via 
the contour integral 

vi(s) = -^i i^^dz, (12) 

where the path of integration, connecting the points —s =p zO on the complex z-plane, 
avoids the cut running along the real negative axis. The integral being traversed in a 
positive (anticlockwise) sense. In this paper we shall assume, without loss of generality, 
that the approximation (|ll|l to the Adler function has only one non-physical singularity 
located on the positive real axis. This is the case, for the exact (explicitly solved) two- 
loop order running coupling in MS like renormalization schemes [j. On the other hand, 
a running coupling at higher orders may be expanded in powers of the exact (explicitly 
solved) two-loop order coupling [61II62| 



(k— loops) /„2n V^ ^(fe) (two— loops)n / „2m /,o\ 

Ots "^ '{Q ) = 2_^Cn al "^ ' (Q ) I exact, (13) 



where the numerical coefficients C„ are determined in terms of the /3-function coef- 
ficients (see Appendix A). It was shown in [33] that this series has a sufficiently large 
radius of convergence in the space of the coupling constants, and its partial sums pro- 
vide very accurate approximations to the exact k-th order (k > 2) coupling in the 
complex Q plane. To construct accurate approximations to the running coupling for 
small values of |(3| , one should keep sufficiently large number of terms in the partial 
sum. The Adler function evaluated with this approximation to the coupling has only 
one non-physical singularity located on the positive Q -axis. The corresponding cut 
runs along the finite interval of the positive Q -axis. Nevertheless, formula (|12p is still 
valid provided that the integration contour avoids the physical as well as non-physical 
cut. 



^ Tlie analytic structure of the explicit exact solution to the RG equation at the two-loop 
order has been determined in |30II31||32) . 




Fig. 1 Contour in the complex Q^ plane used in the Cauchy relation II15II . Branch points on 
the real axis are represented by the blobs and brunch cuts by the zigzagging line. 



Let us separate out the parton level term from the perturbative Adler function 

oo 

^RGl(O') = 1 + dRGl(Q') : dRGl(Q') = Yl dkCisiQ^), (14) 

fc=l 

where as{Q ) — as{Q )/it. As it was discussed above, the function dRGl(Q ) is analytic 
except the cuts running along the real Q -axis. The physical cut runs along the real 
negative semi-axis — oo < Q < 0, and the non-physical cut runs along the positive 
interval < Q < sL) where the point Q = Q^ = ^L > corresponds to the "Landau 
singularity" . We may then write a Cauchy relation 






(15) 



where the integral is taken round the closed contour 7^ drawn in Fig.l. The contour 
consists of the arc of the circle \Q — sl| = sl, straight lines parallel to the real negative 
Q axis and passes round a big circle. Using formula (|15|) together with the asymptotic 
condition dRcK-^^) — >■ as |z| — >■ oo, we derive a violated dispersion relation (DR) 



rfRGl(Q')=dAPT(0')+dL(Q') 



here the function dAPT(Q ) satisfies the DR 

rfAPT(Q^) = : 
with the effective spectral density 



a + Q2 



da. 



Pcsi^) = Im{dRGi(-cr - lO)}. 



(16) 



(17) 



(18) 



It is to be noted here that the function 



^APT(Q') = l + dAPT(Q') 



(19) 



is the analytic image of the perturbative Adler function determined in the sense of the 
Analytic Perturbation Theory (APT) approach of Shirkov and Solovtsov [40II41| . The 
second term in p6p . which violates the DR, corresponds to the contribution to the 
integral (|15[) coming from the "Landau branch cut". It is represented by the contour 
integral 

-,2x _ 1 / rfRGl(C) 



taken round the circle {C • C — ^L + sl exp (z(jf>), — tt < (/> < tt} in the positive (anti- 
clockwise) direction. 

The perturbation theory approximation to the hadronic spectral function is calcu- 
lated by inserting the series ()14|) into the inversion formula H12|) . An important point 
is that the "Landau part" di^{Q ) does not contribute into the spectral function, pro- 
vided that s > 0. To see this, let us evaluate this contribution to the spectral function, 
with the aid of formula H20|l . 



1 /■ ( 1 /-^+'" 1 -I ^ ' 



here we have interchanged the order of integration in the repeated integral. Let us 
consider the integral under braces. For C 7^ the integrand has two simple poles inside 
the contour of integration. It follows from the theorem of residues that this integral 
vanishes, provided s > 0, 

1 r — s+lO ■, 



and the same result holds for ( — 0. We have thus found that only the "analytic 
component" dAPT(Q ) gives a finite contribution into the hadronic spectral function. 
Using DR p7p and inversion formula p2p . one finds the expression for the spectral 
function in terms of the effective spectral density 

.PQC°(.) = .fPT(.) = i(l + r(s)), (22) 

where 

1 f°° Pcff(o') 



r 



(s)^- ^^SEAl^d^. (23) 



Note that formulas (|22p and 1)23^ were previously obtained in the context of APT (see 
[39| and [S]). With the help of formula H22p . we express the "perturbative component" 
of the total "experimental" Adler function in terms of the effective spectral density 



IC{Q^,s){l + r{s))di 



(24) 



where we have introduced the notation JC{Q ,s) = Q /(s + Q ) . Integrating p4p by 
parts we obtain a more convenient representation 
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Let us now evaluate power suppressed corrections to the total "experimental" Adler 
function (|8|). We may rewrite the perturbative component of the Adler function iden- 
tically 

Jo ,^ (26) 

Jo 
in the first line of (|26p we have used the definition of the analytic image of the Adler 
function 

/'OO 

^APT(Q') = l + dAPT(Q') = 2 / IC{Q^s)vt''^{s)ds, (27) 

■Jo 

which is easily deduced from the discussion given above. The last equality on the right of 

(|26p follows from formula p6|) . The power suppressed part of the total "experimental" 

Adler function is determined as 

^p.s.(Q',sp) = ^.'cxp"(0') - Orgi(Q')- (28) 

Combining formulas ([HI, (|16p and (|26p . we rewrite formula (|28|l in the form 

^2 " ^ - Dexp(Q^ Sp) + OpQCD(Q', Sp) - ^RGl(Q') ,„Q^ 

2 „\o„,cxp/„\ . „ I /^2\ fSp r^/^2 „\o„,APT/„\ j„ {^^1 



Dp.s.iQ', Sp) = Dcxp(Q^ Sp) + DpQCD(Q^ Sp) - £)RGi(g 

= /;' i^(Q2, s)2^f P(s)ds - dL(Q2) _ /^''P i^(Q2, s)2<PT(^)^^ 



From definitions 1)20^ and H24|) . we obtain the asymptotic formulas, for Q — >■ oo, 

/C(Q^s)^Q"V0(sQ"'*), di^{Q^)^c-^A^Q-^ + 0{A'^Q-'^), (30) 

where /l denotes the conventional MS-scheme QCD parameter {A = Ayrg). Since the 
parameter sl is proportional to A q the coefficient cl is a positive number independent 
of yl 

CL = ^"'^y'^rfRGl(C)rfC=^^y" dRGl(sL + SLe''^)d<^. (31) 

Using formulas (p9)) and (|30)) . we write asymptotic expansion for Dp.s.{Q ,Sp). It 
follows from the OPE that the leading term in the asymptotic expansion, proportional 
to Q~ , vanishes if the quarks are massless. This leads to the equation [j 

cltI^ + / " 2vt^'^{s)ds = ci^A^ + Sp + I " r{s)ds = / " 2v°''P{s)ds, (32) 

Jo Jo Jo 

the first equality in Eq. ((32]) follows from the relation ((22]). Using Eq. (|23p . by partial 
integration we find 



/■Sp ^ rSj, 

/ r(s)ds = Spr(sp) + - / Peff(o- 
Jo '^ Jo 



)da, (33) 

here we have used the relation sr{s) — > as s — >■ 0, which holds in every order of 
perturbation theory. Combining Eqs. ()32p and (|33l) . we obtain 

CL^l^ + sp(l + r(sp)) + - f " Pes{'y)da = / ' 2vf'^{s)ds. (34) 

^ Jo Jo 



* The expressions for sl in terms of A up to fourth order in perturbation theory may be 
found in I33| (see, also. Appendix A). 
^ The FOPT version of this equation reads Jq '' v^'-' (s)ds = Jq'' f°''^(s)ds (see 3 ). 



11 



Table 1 The numerical values of the coefficient cl in the MS scheme as calculated from 
formula II31I I in the case of the four-loop order /3-function. 





Approximations to the Adler function 




LO NLO N^LO NSLO N^LO 


CL 


0.301262 0.453421 0.555401 0.651373 0.721687 



The coefRcient cl is calculated numerically from formula H31|) . In this calculations we 
use the exact (explicit) two-loop running couphng and exact (numeric) four-loop run- 
ning coupling |_j. In particular, using the next-to-next-to-leading order approximation 
(N LO) to the Adler function constructed in terms of the exact two-loop order running 
coupling we find 

CL!two-loopbcta = 0.421163. (35) 

The numerical values of the coefficient cl evaluated in the MS scheme in the case 
of the four-loop order exact (numeric) running coupling are listed in Table [T] In the 
calculations we have used the approximations to the Adler function of increasing order 

I |. For the unknown 0{ag) correction to the Adler function, we employ the geometric 

estimate d^ = d^^d^^/d^) — 378 [12|. 

It follows from the mixed representation (U) for the spectral function that one may 
calculate in perturbation theory the decay rate of the r lepton into hadrons of invariant 
mass larger than y'Sp 

i?^y''|s>Sp =6lVudPS'EW / wr{s)vt^'^{s)ds, (36) 

■^Sp 

where 

m^ \ rrir J \ m^ 

Vud a-iid Sew denote the flavor mixing matrix element and an electro-weak correction 
term respectively [6]. Equation (|36|l reduces to 



J s„ 



Wt{s)vi {s)ds — I Wt{s)v°^^ {s)d s. (37) 



Using relation (|22p . we express the left hand side of (|37|l in terms of the effective 
spectral density. By integrating by parts, after some algebra, we obtain 






4 



wr{s)vt''^{s)ds = i ( 1 - iiL ) ( 1 + ill ) (1 + risp)) 



mi 



1 nPosis)f, s\ f^^s^^^ 



4:71 L s \ ml I \ m? 



(38) 



To clarify the difference between the APT and APT frameworks, a few comments 
are in order: i) In APT the spectral function is determined through formula (|22p in 



^^ Application of the explicit series solution II13I I for the four-loop coupling yield the same 
results. 

^^ we will use the abbreviation N'^LO to denote the order 0{as ) approximation to the 
Adler function. 
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the entire region < s < oo, whereas in APT this formula holds only for s > Sp. ii) 
In APT the RG improved approximation to the Adler function, -Drgi(Q )i is replaced 
with corresponding "analytic" image -Dapt(Q ) from the outset. Then formula H28p . 
the definition of the power suppressed contributions to the Adler function must be 
suitably modified (see [47]). In this paper, we do not mention this procedure, iii) To 
parametrize our results we use the standard coupling constant as{mr), whereas in 
APT the results are parametrized in terms of the analytic coupling as{m^)a,n 



3 Numerical Results for the Parameters 

To extract the parameters Sp and A from the data we have to solve the system of 
equations 

<Pi{sp,A')= r%fP(s)ds, (39) 

Jo 

/•ml 
*2(Sp, /!')=/ Wr(s)vl''P{s)ds, (40) 

J Sp 

where the functions ^1^2 are defined as 

^l{sp,A^) = ^(1 + r(sp)) + ^ r pM'^)da+ ^A', (41) 

•^ 



$2iSp,A^) = {l-Spf{l + Sp) 



(l + r(sp)) 



4^ / 7. (1-J/) (l + y)rff 

with Sp = Sp/mr- The right hand sides of Eqs. (|39|l - (|40| l are determined in terms of the 
empirical function w™P(s). We reconstruct the experimental vector spectral function 
from the ALEPH 2005 spectral data for the vector invariant mass squared distribution 
which is publicly available [1^ (see Appendix B). The spectral function is measured 
at discrete points of the energy squared. To interpolate the spectral function between 
these points we use cubic splines. 

We solve the system of equations (|39|) - (|40[ ) numerically using various approxima- 
tions to the Adler function. Since the system is transcendental it has more than one 
solution. In Table[2] we present the first reasonable solution for the parameters obtained 
at next-to-next-to- leading order (N LO). From the Table, we see that the predictions 
for Sp are stable with respect to the loop corrections to the /^-function. In this regard, 
the predictions for the QCD scale parameter is more sensitive. The two values of A 
obtained with the two- and four-loop /^-functions differ in about 10%. However, this 
corresponds to the small difference as('7i^)|fojjr-loop ~ '^sirn^)\tvjo-loop ~ 0.0017. 

The solution for Sp obtained with the two-loop running coupling should be com- 
pared with the estimate Sp = 1.60 ± 0.17 extracted in [5] from the earlier ALEPH 
data. Our prediction for the central value of Sp (see Table [2]) is greater in about 
7%. However, with the more accurate data, we have obtained smaller experimen- 
tal errors on the parameters (see Appendix B). Our estimate for the central value, 
^1 {two-loop /3} = 383 MeV is somewhat above the value ^| {two-loop /3} — 372 MeV 
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Table 2 The first solution for the parameters Sp and A = Ajija obtained at N^LO. The two- 
and four-loop running couplings have been used. The extracted values of the strong coupling 
constant as{m^) are also given. The error bars refer to the experimental uncertainty only. 



Observable Approximation to the /3-function 



Two-loop Four-loop 



Sp GeV^ 1.711 ±0.054 1.709 ± 0.054 

A GeV 0.383 ± 0.034 0.348 ± 0.030 

as{m'^) 0.320 ill 0.015 0.321 ili 0.016 



Table 3 The same as in Table [2] for the case of the second solution for the parameters. 



Observable Approximation to the /3-function 



Two-loop Four-loop 



Sp GeV^ 0.606 ± 0.003 0.607 ± 0.003 

A GeV 0.583 ±0.018 0.522 ± 0.016 

asjm'^) 0.417 ±0.010 0.424 ±0.011 



Table 4 Numerical values for the parameters in the MS scheme extracted from the t data 
order- by-order within the modified procedure based on APT+. 



Observable 


Approximation to the Adler function 




LO 


NLO 


N^LO 


N^LO 


N-^LO 


Sp GeV2 
ylGeV 
as (ml) 


1.707 
0.486 
0.401 


1.710 
0.378 
0.337 


1.709 
0.348 
0.321 


1.707 
0.332 
0.313 


1.705 
0.323 
0.308 



accepted in f3l. However, one should keep in mind that in ^Sj only one equation, the 
FOPT counterpart of Eq. H39|l . has been utilized. 

Note that the system (|39p - (|40p permits one more solution for the parameters in the 
range 200 MeV < yl < 600 MeV (see Table El). An attractive feature of this solution is 
that it predicts a smaller value for the onset of perturbation theory: Sp = 0.607 GeV ~ 
rrip {mp stands for the p-meson mass). However, considering current status of as we 
find the extracted value for the strong coupling constant too large. For this reason, we 
decline this solution. 

We also determine the experimental uncertainties on the parameters coming from 
the uncertainties of the vector invariant mass squared distribution. The correlations 
between the errors of the distribution are properly taken into account. Cumbersome 
technical details of the error analysis are relegated into Appendix B. 

It is useful to determine the so-called indicative estimates of the theoretical uncer- 
tainties on the numerical values of the parameters (for the definition see [27]). This 
requires us to test convergence of the numerical results order-by order in perturbation 
theory. We use consecutive approximations to the Adler function from LO to N LO. 
For the unknown 0{as) correction, we use the geometric estimate ds — d4{d4/d3) = 
378 ± 378 [T5]. The results for the extracted values of the parameters are presented 
in Table |4l Formally, we may write a series for the numerical value of the coupling 
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constant as follows 

4 

Qs(m?)lN'iLO = "s("1?)1lO + X] ^k, 

where Z\j, = as(n^r)lNl'LO ^ '3^s('TiT-)|]^k-iLo- Using the numbers listed in Table|3](we 
use abbreviation APT for the modified APT accepted in this paper) we obtain the 
series 

as{ml)\^i^Q = 0.40f - 0.064 - 0.0f6 - 0.009 - 0.005. (43) 

The changes of the leading term induced by the consecutive corrections in the series 
are found to be: 15.9%, 4.0%, 2.2% and 1.2%. It is interesting to compare the series 
H43[) with its counterpart obtained within standard CIPT. Using the standard CIPT 
to analyze the same data (for details see Appendix C) we obtain the series 

asimr)\'^FJ) = 0.485 - 0.095 - 0.023 - 0.013 - 0.007. (44) 

We see that within CIPT the corrections provide slightly larger changes of the leading 
term: 19.6%, 4.7%, 2.7% and 1.4%. One finds that Zifc(CIPT)/Ak(APT+) « 1.2 for 
fc = 1 — 4 . So that the series H43p converges slightly rapidly than the series H44p . The 
indicative estimate of the theoretical uncertainty is determined as a half of the last 
retained term in the series [27j |_J. As pointed out in [27j , the error defined in this way 
is heuristic and indicative. The actual values of the theoretical errors related to the 
uncalculated higher order terms in the perturbation theory series for the decay rate 
might be even larger (see, for example, papers [51l27ll251l53p . In this paper, however, 
we shall consider only the indicative theoretical errors. From the series H43|) . we obtain 
the estimates 

Qs("i?)Inlo = 0.337 ± O.OlGoxp ± 0.032th 
as(m?)lN2Lo = 0.321 ± O.OlGoxp ± 0.008th 
Qs(m?)lN3Lo = 0.313 ± 0.014oxp ± 0.004th 
Qs(m?)JN4Lo = 0.308 ± 0.014cxp ± 0.002th, (45) 

here we have also included the experimental errors. Analogically, from the CIPT series 
1)44 |l . one obtains 

as(m?)|NLO = 0.390 ± O.OUoxp ± 0.048th 
«s(n^?)lN2LO = 0.367 ± 0.009oxp ± 0.012th 
a.(m?)|N3Lo = 0.354 ± 0.008cxp ± 0.007th 
Qs(m?)lN4Lo = 0.347 ± 0.008cxp ± 0.003th, (46) 

The N LO estimates in (|45|) and ()46|) correspond to the central value ^5 — 378. The 
additional theoretical error in the coupling constant induced from the uncertainty in 
the fifth order unknown coefficient {d^ = 378 ± 378) takes the values 0.0045 (~ 1.5%) 
and 0.0065 (~ 1.9%) in the new and standard extraction procedures respectively |_|. 
Comparing the numbers in formulas H45|) and H46|) . we see that the indicative estimates 
of the theoretical error are smaller within the new procedure. In contrast to this, the 
experimental errors on the values of as increases by the factor of 1.76 within the 



^^ In 27 this definition of the uncertainty has been used within FOPT. 

13 With ds = 756, we have obtained Q!s(m?)|^pT+ = 0.3035 and as{m^)\c-[pT = 0.3407. 
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new procedure . It is remarkable that a more reliable estimate of the theoretical error 
presented in [12] is close to the N LO and N LO values of the indicative error given 
in formula H46p . 

Similarly, determining the indicative theoretical errors on the parameter Sp, we find 
the stable results 

SpInlo = 1-710 ± 0.054exp ± 0.002th GeV^ 

SpIn^lo = 1-709 ± 0.054exp ± 0.001th GeV^ 

SpIn^lo = 1-707 ± 0.054cxp ± 0.001th GeV^ 

spIn-ilo = 1-705 ±0.054cxp± 0.001th GeV^. (47) 

Notice that the ratio as(sp)/as(m.r) ~ 1-22 is not large. However, the APT ex- 
pansion formally depends on the small energy scale v/sp ~ l.SlGeV. So, it is reasonable 
to justify the applicability of the perturbation theory in the APT framework. The 
issue of the applicability of perturbation theory in t decays has been previously ad- 
dressed in 63 . It was pointed out 163, that this question is phenomenological one, and 
it cannot be answered yet from theoretical grounds. In particular, the decay rate of the 
r lepton into hadrons of invariant mass squared smaller than sq (so < it^-t) ha-s been an- 
alyzed within FOPT. Using the ALEPH spectral data, the authors of [63] have deduced 
that the rate can be calculated in pQCD with high accuracy for sq > Smin. ~ 0.7 GeV . 
Note that our estimate for Sp clearly satisfies this condition, Sp/smiu ~ 2.4. Never- 
theless, it is desirable to investigate numerically the convergence of the perturbative 
expansion within APT . Let us derive the expansion for the r- lepton decay rate from 
formula p6|) . The integral on the right of ()36p can be approximated by a non-power 
series. To derive the non-power series, we express the spectral function in terms of the 
effective spectral density using formulas (|22p and H23p . Then we expand the function 
pGff(s) in perturbation theory using formulas H14|) and H18|) . So, we obtain 

5 

^P V*'U>-P = <,"* ■U>-p/{6iV;d|SEw} = E rffe2tfe(m2, Sp) (48) 



where 



k=0 



2to(m?,Sp) = /(sp/m?), (49) 



2 \ / NT/ / 2 



^k>l(mr,Sp) ^ rk(sp)f{sp/mr) / —pk[(T)d(j, (50) 



TT 



2> 



1 / " f{<^/mr 



here we have used the notations; f{x) = -^(1 — x) [1 + x) and 



Pfc(o") = Im{as (-0- - lO)}, 



1 /""" pfc(g-) 

■^ Js 



r-kisp) = - / '^^^da. (51) 



The first term in the series (|48l) . 2lo, corresponds to the (modified) parton level contri- 
bution to the rate. We calculate the functions 21^ numerically by using analytic expres- 
sions for the functions Pfc(o") (see formula (|63p in Appendix A). In the calculation, we 
employ the four-loop running coupling. For the parameters Sp and A = Ajj^, we use 
the numerical values from the Table [IJ namely, the values extracted from the ALEPH 
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Table 5 Comparison of the expansion functions 2t(;(mJ,Sp), A-ki'fTi'^) and the powers of the 
"couplant" as{m^). The four-loop "couplant" is calculated using the value A = 0.3225 GeV. 
To calculate the functions 2lj.(m^, Sp), we have used the values ^|iv[4lo ~ 0.3225 GeV and 
Sp|i^4Lo ~ 1.7053 GeV obtained within APT"*". To calculate the functions y4fe(m^) we have 



used the value A\ 



N-^LO 



: 0.395 GeV obtained from the ALEPH data within CIPT. 



k 


«JH?) 




A-(m?) 


2tfc{m? 


Sp) 


1 


0.9797 


10-1 


0.1511 


0.3275 


10-2 


2 


0.9599 


10-^ 


0.1876 10-1 


0.2400 


10-3 


3 


0.9405 


10-3 


0.2000 ■ 10-2 


0.1455 


10-4 


4 


0.9214 


10-4 


0.1834-10-3 


0.6733 


10-'^ 


5 


0.9028 


10-^ 


0.1383-10-'' 


0.1599 


10-^ 



data within APT at N LO. Using analytically known coefficients dj., k — 
the estimate d^ = 378, we obtain from Eq. H48|l the expansion 



4 and 



pport. 



s>s^ = 0.3747 
+ 0.3304 - lO"'' 



- 10-1 + 0.3275 - 10-2 + 0.3937 ■ 
+ (0.6047- 10"*^) ^0.04127. 



10" 



0.9270 ■ 10" 



(52) 



Consider now the non-power expansion for the perturbation theory correction 5^ ' 
obtained within CIPT [24j (see Appendix C) 



fc=l 



(53) 



where 



Akirnr) 






Re{(l 



l(p 



){l + e"''fa';{mU"''')}d(l> 



to calculate these functions numerically, we employ for the scale parameter A the 
numerical value extracted from the ALEPH data within CIPT at N LO (see Table 
1141) . At N LO, the expansion H53p can be rewritten as 



(5[?j^ =0.1513 + 0.3081 -10"^ 



0.9012 - 10^2 + (0.5233 -10 ^) ^ 0.2091. 



0.1276-10" 

(54) 

Comparing the numerical expansions in Eqs. H52|l and (|54ll . one sees that the APT 
series (I52p displays a faster convergence. In the CIPT expansion (|54p . the corrections 
provide a 38% change of the leading term. In contrast, in the APT expansion H52p the 
corrections provide only a 16% change of the leading term (we recall that the leading 
QCD correction in 1)52^ is the second term in the series). 

The rapid convergence of the series (|52l) may be explained due to the specific 
properties of the expansion functions 2lj,(m,-, Sp). The set of functions {2lj.(m^, Sp)} 
can be viewed as a generalization of the analogical set of functions considered in the 
Shirkov-Solovtsov APT (for properties of the APT expansion functions see [42])- In 
Table [SI we have compared functions 2ti.(m^, Sp) with the functions Ak{rn^). For the 
sake of comparison, we also include in the Table the powers of the "couplant" as('7i,-). 
It is seen from the Table that the functions 2lfe(m,-,Sp) decrease with k much more 
rapidly than the functions Akimr) and ag{mT). 

Usually, it is convenient to perform evolution of the as results to the reference 
scale Mz = 91.187 GeV. This is done by using RG equation and appropriate matching 
conditions at the heavy quark (charm and bottom) thresholds (see [64] and literature 
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Table 6 Estimates for as{M^) obtained from the ALEPH t lepton decay vector data order- by- 
order in perturbation theory. The results obtained within APT+ and GIFT are compared. Two 
errors are given, the experimental (first number) and the error from the evolution procedure 
(second number). 

Perturbative order as(Mf )|^prp+ as(M| )|cipt 

N^LO 0.1187 ±0.0019 ±0.0005 0.1238 ± 0.0009 ± 0.0005 

N^LO 0.1176 ±0.0018 ±0.0005 0.1224 ± 0.0009 ± 0.0005 

N*LO 0.1170 ±0.0018 ±0.0005 0.1217 ± 0.0009 ± 0.0005 



Table 7 Comparison of the RSI and APT+ determinations of the MS coupling constant from 
the T-decay data. Experimental errors are given only. 



Perturbative order 


cs{m^r)\r^^ 


Mn^Dir^" 


NLO 

N^LO 

N^LO 


0.278 ± 0.003 
0.319 ±0.004 
0.312 ±0.004 


0.335 ±0.016 
0.321 ±0.016 
0.313 ±0.014 



therein). The three-loop level matching conditions in the MS scheme were derived in 
[65| . In this paper, we follow the work 66 , where a very accurate analytic approx- 
imation to the four-loop running coupling was suggested. We perform the matching 
at the matching scale m^i^ = 2/iij where /xjj is a scale invariant MS mass of the 
heavy quark ^^i ~ "^h(Mh)- We assume for the scale invariant MS masses the val- 
ues /ic = l-27l:o;JJ GeV and pb = 4.20^0.^7 GeV [67]. Following [66], we evaluate the 
central value and error of as{Mz) according to the formulas 

as{M^) = (a+(Af2) + a7{M^))/2 and Aas{M^) = {at{M^) ~ olJ{mI))/2 

where afiMz) denote the values obtained from af{rn^) = oisirn^) ± Aas{my). In the 

evolution procedure, we have used the exact numeric four-loop running coupling] ]. In 

Table [6] we compare the estimates for as{Mz) obtained within the new (APT ) and 
standard (CIPT) procedures. 

Finally, for the sake of comparison, let us extract the numerical values for the cou- 
pling constant from r decay (V+A) data using the renormalization scheme invariant 
extraction procedure (RSI) of [27]. This procedure is based on FOPT. For the exper- 
imental value of the perturbative part of the r decay rate in the non-strange (V+A) 
channel, we assume the updated value, presented in |26| . 

5'£^p\v+A = 0.2042 ± 0.0050cxp. 

For consistency reasons we use the MS scheme /3-function to the k-loop order with 
the N ^ LO approximation to the Adler function. In Table ([7]) we compare numerical 
values for OLs{mr) obtained within the two approaches, RSI and APT . The relevant 
channels which have been used to extract the coupling are indicated by subscripts. 
One sees from the Table, that beyond NLO there is a good agreement between the two 
methods of the 05(771^-) determination. 



^* We have confirmed that the approximate analytical coupling derived in 1661 leads practi- 
cally to the same numerical results. 
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Table 8 Comparison of the "experimental" Adler function -D"exp" (Q^) with its QCD compo- 
nent DpQcuiQ^ ) Sp) at low momenta. The perturbative component is evaluated within APT+ 
at N^LO using the four-loop running coupling. The absolute and relative statistical errors of 
the "experimental" Adler function are tabulated. 



QGeV 


i?"oxp-(Q") 


OpQCD(Q^Sp) 


't(£'"cxp") 


rel.err. 


0.1 


0.0649 


0.0063 


0.0061 


9.5% 


0.2 


0.2300 


0.0249 


0.0198 


8.6% 


0.3 


0.4354 


0.0545 


0.0333 


7.7% 


0.4 


0.6320 


0.0933 


0.0426 


6.7% 


0.5 


0.7944 


0.1391 


0.0473 


6.0% 


0.6 


0.9162 


0.1895 


0.0484 


5.3% 


0.7 


1.0016 


0.2426 


0.0471 


4.7% 


0.8 


1.0583 


0.2965 


0.0445 


4.2% 


0.9 


1.0942 


0.3497 


0.0412 


3.8% 


1.0 


1.1157 


0.4013 


0.0377 


3.4% 



4 Numerical Results for the "Experimental" Adler Function 

Looking at the numbers in Table [5] we see that our estimates for the parameters are 
somewhat different than those used previously in ^^ . Hence, it is sensible to recalculate 
the experimental Adler function in the infrared region. Another reason to do this is 
the appearance of the improved r data [T^ . More importantly, it is desirable to carry 
out the error analysis too. Furthermore, in contrast to [3], in our calculations we will 
employ APT+. 

The "experimental" Adler function and its QCD component are tabulated in Ta- 
ble |8l The QCD component of the "experimental" Adler function is calculated numer- 
ically at N LO from formula ([25]). In the calculations we employ the four- loop running 
coupling. For the parameters Sp and A, we use the values from Table [5] The absolute 
(±l(j) and relative (in percents) experimental errors of the "experimental" Adler func- 
tion are also tabulated. The error analysis is described in Appendix B. We see from 
the Table that the pQCD component has sizeable contribution to the total "experi- 
mental" Adler function. This contribution increases monotonically with Q from 10% 
(at Q = 0.1 GeV) to 36% (at Q = 1 GeV). 

To test the stability of the numerical results with regards to the higher order cor- 
rections to the /3-function, we have compared two results for the "experimental" Adler 
function that are obtained with the two- and four-loop exact running couplings. The 
pQCD component of the Adler function has been evaluated within APT^ at N LO. 
For the parameters Sp and A, we have used the central values given in Table [21 In the 
region Q = — 1.5 GeV, the difference between using the two- or four-loop approx- 
imation to the /3-function is found to be quite small (~ 0.05%). The approximation 
corresponding to the two-loop running coupling takes slightly large values. 

To test the stability of numerical results with regards to higher order corrections 
to the Adler function, we use various approximations to the pQCD component (see 
Table [5]). We see from the Table that the differences between the consecutive approx- 
imations to the "experimental" Adler function slowly increase as a function of the 
scale. Already, the leading order approximation provides a very accurate result. At 
Q = 1.5 GeV (where the changes induced by the loop correction take maximal values) 
the differences between consecutive approximations (i.e. the differences between the 
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Table 9 Different approximations to the "experimental" Adler function as a function of 
the scale. The function D\ n(Q^) has the pQCD component evaluated within APT"*" at 
NC'^^JLO. To construct this component, we employ the four-loop order running coupling. 



Q GeV 



cxp 



(Q^) 



D?^ „(Q2) 

"cxp ''^ ' 



"exp ' ^^' > 



dI*) 



.(Q^) 



dP „(Q2) 

cxp "^^ -' 



0.1 


0.06494 


0.06494 


0.06494 


0.06494 


0.06494 


0.2 


0.23003 


0.23005 


0.23004 


0.23003 


0.23002 


0.3 


0.43541 


0.43546 


0.43544 


0.43541 


0.43539 


0.4 


0.63196 


0.63205 


0.63201 


0.63196 


0.63192 


0.5 


0.79431 


0.79444 


0.79438 


0.79430 


0.79424 


0.6 


0.91613 


0.91631 


0.91623 


0.91612 


0.91604 


0.7 


1.0015 


1.0017 


1.0016 


1.0015 


1.0014 


0.8 


1.0582 


1.0585 


1.0583 


1.0582 


1.0580 


0.9 


1.0940 


1.0943 


1.0942 


1.0940 


1.0938 


1.0 


1.1154 


1.1158 


1.1157 


1.1154 


1.1152 


1.5 


1.1321 


1.1327 


1.1324 


1.1320 


1.1317 



N^^^^LO and N'^LO approximations) take the values 0.05%, 0.03%, 0.03% and 0.03% 
for fc = 1, 2, 3 and 4 respectively. 

It is instructive to investigate numerically the convergence property of the non- 
power series for the perturbation theory component -DpQCD(Q i^v)- The non-power 
series is obtained from formula (|25p by using perturbation theory expansions for the 
function /9off(s) and r(^s) (see formulas p4|l . ([TS} and (|23p '). The non-power series read 



■DpQCD(Q^,Sp) = ^dfeDfe((5^,Sp), 
fc=0 



(55) 



where 



So(0^ 



i*p; 



(Sp + Q2)' 



Sfe>i(Q ,sp) = 



(Sp + Q2 



■'^fe(sp) 







31 r Pk{<y) 



-da. 



(56) 

(57) 



the functions pk(p) and r;;(sp) are defined in Eqs. (|5ip . Let us truncate the non-power 
expansion (|55|l at N LO (i.e. for fe = 5). Using the N LO estimates for the parameters 
given in Table |3J we evaluate the ratios of the consecutive terms of the series 

7^fc(Q') = (dfe/dfe_l)^fc(Q^sp)/2)fc_l(o^sp), 

for fc = 1 — 5. In Table [TOl we tabulate numerical values of these ratios in the region 
Q = 0.1 — 1.5 GeV. It is seen from the Table, that the magnitudes of the ratios are 
sufficiently small to guarantee fast numerical convergence of the series: TZ^iQ ) ^ 0.368 
(fc = 1 — 5) for all values of Q in the considered interval. 

Let us now compare our numerical results on the "experimental" Adler function 
with the previous results of work ^3 . First, we repeat the calculation within the ap- 
proach of [3] using the improved data. Assuming the value ^ttf- = 372 ± 76 MeV used 
in [3], we solve numerically the FOPT counterpart of the equation H39|l . Thus we find 
the solution Sp = 1.621 ± 0.163 GeV . The central value of this estimate is slightly 
large, by 0.021, than the value obtained in [3]. Using the values Sp = 1.621 GeV and 
^MS ~ ^^^ MeV, we calculate the "experimental" Adler function within the (modified) 
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Table 10 The ratios of the consecutive terms in the non-power series II55II as a function of 
the scale. 

QGeV 7^l(Q2) ■RiiQ'^) Tl^iQ'^) UiiQ'^) TZsiQ^) 



0.1 


0.077 


0.112 


0.228 


0.366 


0.238 


0.3 


0.077 


0.111 


0.228 


0.366 


0.240 


0.5 


0.077 


0.111 


0.227 


0.366 


0.242 


0.7 


0.076 


0.110 


0.226 


0.367 


0.245 


1.0 


0.075 


0.109 


0.225 


0.367 


0.251 


1.3 


0.073 


0.108 


0.223 


0.368 


0.257 


1.5 


0.072 


0.107 


0.222 


0.368 


0.261 



Table 11 Comparison of the "experimental" Adler functions evaluated within the modified 
FOPT and APT+. The pQCD components of the functions are constructed at N^LO with 
the two-loop order running coupling. The pQCD component of D{Q'^)\f,^^'^ corresponds to 
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the values ylj^ = 372 MeV and Sp = 1.621 GeV^. The pQCD component of D{Q^)\ff^'-„ 

corresponds to the values Ajj^ = 383 MeV and Sp = 1.711 GeV^. The relative difference 
between these functions is also tabulated. 



QGeV 


D{Q')\^S.7 


D(Q2)|APT+ 


rel.diff. 


0.1 


0.0651 


0.0649 


0.31% 


0.2 


0.2305 


0.2301 


0.20% 


0.3 


0.4364 


0.4355 


0.21% 


0.4 


0.6336 


0.6321 


0.24% 


0.5 


0.7967 


0.7945 


0.28% 


0.6 


0.9191 


0.9164 


0.29% 


0.7 


1.0050 


1.0018 


0.32% 


0.8 


1.0621 


1.0586 


0.33% 


0.9 


1.0982 


1.0945 


0.34% 


1.0 


1.1197 


1.1160 


0.33% 


1.1 


1.1316 


1.1280 


0.32% 


1.2 


1.1371 


1.1337 


0.30% 


1.3 


1.1386 


1.1355 


0.27% 


1.4 


1.1377 


1.1350 


0.24% 


1.5 


1.1354 


1.1330 


0.21% 



FOPT at N LO. This should be compared with the new approximation computed in 
the same order within APT . In the case of APT , we use the values A-rj^ — 383 MeV 
and Sp = 1.711 GeV (see Table [2]). To be consistent with [2, we use the two-loop 
exact running coupling. In Table 1111 we compare numerically two approximations to 
the "experimental" Adler function, the functions £)™p,7'(Q^) and -D^xp" iQ^)- From 

the Table, we see that the functions are close, but Du^ ^^ {Q ) > D,,^ „ {Q ). The 
relative difference between the functions in the considered region varies in the interval 
0.20% - 0.34%. 



5 Conclusion 

We have extracted the numerical values of the strong coupling constant as and the 
parameter Sp (the square of the boundary energy) from the non-strange vector r data 
provided by ALEPH. Based on the semi-empirical representation Q for the hadronic 
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non-strange vector spectral function, we have developed a modified extraction pro- 
cedure. This procedure enabled us to avoid direct application of the standard OPE 
formalism in Minkowski space. The distinguishing feature of our analysis is that we 
have determined the two parameters {as and Sp) simultaneously from the data. 

In Sect. 2, we have derived a violated DR for the RG improved perturbation theory 
correction to the Adler function, the formula H16|) . Using the violated DR, we have 
shown that the perturbation theory component of the "experimental" spectral function 
is determined via the APT formula H22p . This determines the hadronic spectral function 
in terms of the effective spectral density, poff (''")i the basic object of the perturbation 
theory calculation. We have obtained a convenient expression for the pQCD part of the 
"experimental" Adler function in terms of the effective spectral density, the formula 
(|25|) . Using the violated DR (|16l) . we have determined the power suppressed corrections 
to the "experimental" Adler function via the formula H29|l . Making further use of the 
consistency condition from the OPE for the "experimental" Adler function, we have 
derived Eq. (|34p . This equation relates the parameters Sp and A to the values of the 
hadronic spectral function on the range < s < Sp. Next we used the ansatz Q for the 
spectral function to calculate the r decay rate Rt,v\s>Sp- In this way, we have derived 
Eq. (|37|) which relates the parameters to the integral of the hadronic spectral function 
(multiplied by known function) over the range Sp < s < 111^. 

In Sect. 3, we have solved, numerically, the obtained system of equations for the 
parameters Sp and A = Aj-^. To examine the convergence of the numerical results 
for the parameters, we have used perturbation theory approximations to the Adler 
function up to the N LO. The indicative estimates of the theoretical errors [27] are 
used as a criterion of the quality of the approximations. Based on this criterion, we have 
demonstrated that the new framework (APT ), compared to the standard one (CIPT), 
provides a better numerical convergence for the extracted value of the coupling constant 
as(TO-,-). It is remarkable that the central values of the coupling constant extracted 
within APT in different orders of perturbation theory become systematically smaller 
as compared to the corresponding values obtained within CIPT (cf. formulas (|45p and 
(|46|l '). The changes in the central values are not within the quoted experimental and 
theoretical errors. At N LO The central values of asimT-) in formulas (I45p and (|46[) 
differ from each other in about 2.7 standard deviation, if the error is determined within 

APT+, a = i/o-lxp. -I- o-2j^ ^ 0.0151. With the error obtained in CIPT, a « 0.0107, 

one finds even large difference, 3.8 a I I. 

We have examined the stability of the numerical value for the duality point s^ with 
regard to perturbation theory corrections. We have obtained a surprisingly stable result 
(see formulas (|47|) ) s^ = 1.71±0.05cxp . ± 0.00th. GeV . Our prediction for the central 
value of this parameter is higher in about 7% than the value presented previously in 

m- 

Having included into analysis the fourth order coefficient d4, we achieved excellent 
agreement between the lattice and tau-decay determinations of the strong coupling 
constant (at N LO the central value of the constant given in Table (O coincides with 
the central value quoted in 0). For this reason we believe that APT^ provides better 
approximation than CIPT. 

To justify the applicability of APT in calculations of the r-lepton decay rates, 
we examine numerically the APT"*" series. The APT+ expansion for the rate Rt.v\s>Sp 
represents asymptotic expansion over a non-power set of specific functions {2ln (mr , Sp ) } 



Due to the larger experimental error obtained within APT+, o"apt/o"CIPT ~ 1-4. 
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rather than the powers of as{mT). The APT and CIPT series for the rates -Rr,v|s>Sp 
and Rt V have been compared numericaUy. We have confirmed that the APT"^ expan- 
sion displays a faster convergence. 

Our approach confirms that there is a theoretical systematic uncertainty not in- 
cluded in the error assessments obtained in previous studies by ignoring the higher order 
OPE contributions, the conclusion achieved in work [ij. In this connection, our study 
suggests that the truncated OPE series cannot approximate sufficiently accurately the 
integrals of the spectral function over the low energy region < s < Sp ~ 1.7 GeV . 

In Sect. 4, we have recalculated the "experimental" Adler function within the APT^ 
prescription using the new estimates for the parameters Sp and A. In addition, we have 
determined the errors on this function coming from the uncertainties of the parameters 
and spectral function. Numerical results for the Adler function obtained within APT 
have been found to be remarkable stable in perturbation theory (see Tables [9l and [TO] 

)■ 

In Appendix A, we have given practical formulas for numerical calculation of the 
MS running coupling at higher orders. The Lambert-W solutions to the RG equation is 
reviewed. An accurate analytic approximation to the effective spectral density Pes{cr) 
at higher orders is derived. In Appendix B, we have derived formulas (within APT ) 
for calculating the experimental uncertainties on the extracted values of the param- 
eters and on the "experimental" Adler function. In Appendix C, we have analyzed 
the ALEPH non-strange vector spectral data within the standard CIPT prescription. 
Namely, we have performed some necessary calculations needed for comparing the 
CIPT and APT+ prescriptions (see Sect.E)). 

The procedure suggested here can be obviously extended to analyze the non-strange 
r-data from the axial-vector (A) and vector plus axial-vector (V-l-A) channels. To check 
the reliability of the new extraction procedure, it is desirable to compare the V, A, and 
V-l-A determinations of the coupling constant. A similar framework may be constructed 
on the basis of FOPT. This will enable us to estimate total theoretical errors on the 
extracted values of the parameters A and Sp. It should be remarked that a shortcoming 
of the ansatz Q is that it completely ignores the non-perturbative contributions to the 
spectral function coming from the region s > Sp. The importance of these contributions 
for accurate determination of the coupling constant has been demonstrated in recent 
studies |20II21) . We hope to report on these aspects in future publications. 
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A A Series Solution to the Renormalization Group Equation 

In our notation the RG equation for the running coupling reads 

^^a, =/3(a,) = -a2^/3„aJ, (58) 

ri — 

where as = as{Q'^) = Os(Q^)/7r with Os(Q^) being the running coupling. In the MS scheme, 
the /3-function coefficients are known to four loops '69'. For three active quark flavours the 
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first four coefficients take tlie values 

/3o = 9/4, /3i =4, 132 = 10.05990, ft = 47.22804. 

In general, the RG equation II58I I. to an arbitrary order in perturbation theory, can not be solved 

explicitly for the coupling. Usually, the equation is solved in the asymptotical regime ^p" S> 1. 
For our purposes the asymptotic solution is not suitable, since we need an accurate solution at 
relatively low energies. One may, of course, solve the RG equation numerically. However, it is 
more convenient to derive some accurate analytic approximation to the coupling. Fortunately, 
it is possible to solve the RG equation explicitly for the coupling at the two-loop order 1301 
I31| . The explicit expression for the MS scheme running coupling at the two-loop order reads 

where /3o and /3i are the first two /3-function coefficients 

1 / 2 \ 1 / 38 

/3o = - (ll - -n,) , ;3, = -(l02»-n, 

bl = /3i//3g, A = Ajja and W—i denotes the branch of the Lambert W function [70J. 

The coupling to higher orders may be expanded in powers of the two-loop order coupling 

mm 

ai'=>^)(Q2) = ^,W42)"(Q2), (60) 

n = l 

The first two coefficients in this series are universal: c^ = 1 and Cj = (the condition 

(k) 

Cj = follows from the conventional definition of the A parameter). Other coefficients are 
determined in terms of the /3- function coefficients. The four-loop expressions for the first several 
coefficients are given by 

(4) _ & (4) _ ft (4) _ 5 /fty ftft 



^ ft' * 2ft' 5 3 Vft/ 6/32 

It was proved in |33l that the series has a finite radius of convergence, and the radius is suf- 
ficiently large for all Uf values of practical interest. Partial sums of the series 16011 provide 
very accurate approximations to the higher order coupling in the wide range of Q^. In par- 
ticular, these approximations may be safely used at low energies. Thus, for Q = 1 GeV and 
A = 0.347 GeV, the partial sum with the first twelve terms reproduce the exact four-loop 
coupling with the precision better than 0.02%. Using the exact solution II59II . the analytical 
structure of the two- loop coupling in the complex Q^.pi^ne has been determined I30II31| . It 
was found that the coupling is an analytic function in the whole complex plane except the cuts 
running along the real Q^ axis. Besides the physical cut {Q^ : —00 < Q^ < 0} corresponding 
to the logarithmic singularity at Q^ = 0, there is also the "Landau" cut {Q^ : < Q^ < Qj^} 
corresponding to the Landau singularity on the positive Q^.^xis. The Landau singularity is 
a second order algebraic branch point located at Qj^ = b^ ^ ^ (^r ^ ~ 1.205 for Uf = 3). 
The relevant branch of the Lambert function on the complex Q^ plane is determined by the 
analytical continuation. For the physical vales of Uf {0 < Uf < 6) the relevant branch is W—i 
on the upper-half plane, whereas the branch is Wi on the lower-half plane. A limiting value of 
the coupling from above the physical cut (Q^ = _o- -|- lO, c > 0) is then determined by I32| 

al i—u + tU) = ; r with C+ = -— exp < — jtt 1 

^ ' ftl-fiy_l(C+) ebi\A^J "^X \bi 



similarly, one may write 



(61) 



4^\-a-^0) = -^ ^—— with C- = — (^Y'^''e.p\^J^-l]}. (62) 
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Table 12 Numerical values of the first twelve coefficients in the expansion II6OII for the four- 
loop running coupling at nj = 3 quark flavours. 



n 


c„ 


1 " 


Cn 


1 


1 


7 


392.1241 


2 





8 


2413.463 


3 


3863/864 


9 


8248.857 


4 


10.49512 


10 


31348.18 


5 


27.09804 


11 


147697.8 


6 


190.2642 


12 


507565.0 



Note that the limiting values a^ (—a it «0) satisfy the Schwarz 'principle of reflection', 
tts (— cr — lO) = cVa (— o" + *0), provided that W has near conjugate symmetry I70| : 

Wfc(^) = VK_fc(2). 

We may construct an analytic approximation to the Adler function in perturbation theory 
using formula II6OII for the four-loop running coupling. In this approximation, the Adler function 
is an analytic function in the cut complex Q^ plane. It has the branch points at Q^ = and 
Q^ = Qj^ = ^i ^ A^ > 0. The effective spectral density III8I I associated with the Adler function 
is then readily calculated, leading to the analytic expression 

PeffH = Im|^d„( f^c<;t'4'''"(-a-i0)) I, (63) 

where N > 3 and a), {—cr — lO) is determined in terms of the W function as given in formula 
II62II . Formula II63II considerably simplifies numerical calculations of integrals of the effective 
spectral function (computer algebra system Maple has an arbitrary precision implementation 
of all branches of the Lambert function). In the most of the calculations, we have used the 
truncated series II6OII for the four-loop running coupling preserving the first twelve terms in the 
series. Numerical values of the first twelve coefficients of the series, for Uf = 3 quark flavours, 
are tabulated in Table [121 



B The Error Analysis 

In this appendix we will evaluate the experimental errors on the extracted values of the pa- 
rameters. We will also determine the errors on the "experimental" Adler function. The main 
quantity employed in our analysis is the vector (non-strange) spectral function vi{s). It is 
related with the vector invariant mass squared distribution (the function sfm2{s) in the no- 
tations of |10p 

vi(s) = K{s)sfm2(s), (64) 

the kinematical factor re(s) is 

I \ — \r "^? f ^^\ ^ 

'^^'' ~ (6|V;dp5Ew) V bT; (l-s/m2)2(H-2s/m2)' ^ ' 

where [Vmj| = 0.9746 it 0.0006 denotes the flavor mixing matrix elements, the factor Sew = 
1.0198 it 0.0006 is an electro-weak correction term, nir = 1777.03l[]'^g MeV, By = (31.82 it 
0.22)% and Be = (17.810it0.039)% are the vector and leptonic branching fractions respectively 
(in this paper, we assume these estimates following |9]), A/" is the normalizing constant 



Af = I " sfm2{s)ds \ 



1 



0.794748 



The quantity sfm2{s) is measured at 140 equidistant values of the energy squared variable 
starting from si = 0.0125 GeV with the bin size Z\i,i„ = 0.025 GeV^. Note that the factor k(s) 
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Table 13 A few measured values of sfm2{s) and vi{s). The standard errors for these quan- 
tities are indicated. 



sGey2 


sfm2{s) 


0"stm2(s) 


Dl(s) 


0-1,1 (s) 


0.0875 


0.004923 


0.001251 


0.006027 


0.001531 


0.1125 


0.022630 


0.003092 


0.027745 


0.003791 


0.1375 


0.037048 


0.004520 


0.045504 


0.005551 


0.1625 


0.056542 


0.005747 


0.069597 


0.007074 


0.1875 


0.073407 


0.005875 


0.090583 


0.007250 


0.2125 


0.095429 


0.006541 


0.118095 


0.008095 


0.2375 


0.122440 


0.007574 


0.152005 


0.009403 



is determined within an accuracy of less than 1% for all values of s in the range s = — m^, 
while the errors in determination of s/m2(s) are considerably large. Hence, it is safe to ignore 
the uncertainties coming from the factor k{s). We may then write 



0-1,1 W = \l^{^k)\o-sfm2[k], fc = 1, ...140 



(66) 



where o^^ [A;] and cr^j,-^2lf^] stand for the standard deviations of 1)1(84.) and sfm2{sk) respec- 
tively, and Sfc = si + (fc — l)Ai,i„ (k=l,2. . . ). By definition 



at^[k]=E[{^i(st:)-vi{s^)r], 



(67) 



etcF^. Similarly, the covariance matrices of the errors associated with the quantities vi{s) and 
sfm2{sf;) are related by the formula 



Cik = K{Si)K(s^)Cik, 



(68) 



where Cik = cov{vi{si),vi{si^)) and C^fc = cov(sfm2(si), sfm2{si^)). So that respective cor- 
relation coefficients coincide 



o-„l[A:]o-„i[/] crsfm,2lk]o-sf,„2ll] 



(69) 



In Table [131 we present a few measured values of sfm2{s) and vi{s) together with the as- 
sociated uncertainties. Our goal is to estimate the uncertainties on the extracted values of 
the parameters induced from the experimental uncertainties of the spectral function. We start 
from the system of II39I I-I I40I1 . which we rewrite in the form 

<Piix,y)=£i{x,{vi}) (70) 

■P2{x,y) = £2{x,{vi}) (71) 

where we have introduced the notations x = Sp, y = A^ and 

ei{x,{vi}) = f vi{t)dt 

Jo 

[■ml 
£2{X,{V^}) = / Wr{t)vi(t)dt, 

J X 

to avoid a cumbersome notation the superscript "exp." in function «^''^(s) has been omitted. 
The solution to the system II70II - II71II should be considered as a functional of vi{s). Let a 
solution for the parameters, for a given function vi{s), is 



X = i/JidiJi}) 
y = 'i/'2({«l})- 



(72) 
(73) 



The symbol E[x] refers to the mean value of x. 
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we may write vi{s) = vi{s) + Svi(s) , where vi{s) is the central (average) value and Svi(s) 
is the deviation. The central values of the parameters should be determined by solving the 
system Il70ll - (l71l l for vi{x) = vi{x) (see, for example, the book 1681 ) i.e. 

x = iPi{{vi}) (74) 

y = M{vi}). (75) 

Let us expand the functionals £1 2(2;, {tii}) in powers of a small variation 5vi(s), preserving 
the terms linear in Sx and 5vi{s) 



Si 
£2 



(x,{vi}) = <^i{x,y) + Sxvi{x)+ I 5vi(t)dt + ... (76) 

Jo 

(x, {Vl}) = ^2 {x,y) - W-r {x)vi (x)Sx + Wr (t)5vi (t)d t + ..., (77) 



here use has been made of the equations £1^2(2;, {1)1}) = 'Pi^2{x,y). Insert these expansions 
into Eqs.l l70l l- (I71| | and expand the left hand sides of the equations in powers of 5x and 5y 

d<Pi2(x,v)r. d^i2(x,y)^ 
■^iM^, y) = -Piais:, y) + 'll '"' 5x + ^'f} ''^' Sy + .... 

ox ay 

Retaining terms linear in 5x, 5y and Svi, we are led to the following linear algebraic system 
of equations for the variations Sx and Sy 



(78) 
where 





AiSx + BiSy = Gi 




A2&x + B2Sy = G2, 


_ d'Pi(x,y) 
dx 


-.,(.), Bi-'^f'y^ 

dy 


'^-^''y^ ,.,f. 


V.,.., r,^_d^2(x,y) 



■■ / Svi{t)dt 
Jo 



dx dy 

Using the explicit formulas 1141 l l and I I42I I. after some algebra, we obtain 

d'P]_{x,y) _ {l + r(x)) 

dx ~ 2 

d^l{x,y) 1 /-inC^/a) 



G2= Wr(t)Svi{t)dt. 

J X 



dy 



I rl-n{x/y) 

— p,s{t)e'dt+-^ 

2vr 7_oo 2 



d'p2{x,y) _ (l + r(a:))^ f x 



dx 4m,| V 'TtJ 

/'ln(i 

^ / Pcff(i)e'P 

^T Jln{x/ 



d<P2{x,y) ^ 1 /■'"''"-''''' -^^(jyp/?^ 
dy 4nml Jin(x/y) ° V'"!? 



where Peff (*) = Poff ("") with cr = A^ exp(i), P(z) = 2(z — 1)^(22 + 1), and r{s) is defined in |[23ll . 
After solving the system jTSj , we take the averages of the deviations squared (the variances) 



(Sxf = (BlGl + BlGl - 2BiB2GiG2)/V'^ 

(5yf = {AiGl + A^Gf~2AiA2a[G^)/V^, (79) 



where T> = A1B2 — A2B1 , and the overlined symbols refer to the averages: (Sx) = E[{x — x)^] 
etc. To calculate the averages G|, G| and G1G2 we replace the integrals Gi,2 by sums over 
the equidistant mesh, using the trapezoidal rule, 

Gi^Aj29kSvi{tk), G2 ~ A J2 ^kSviit^:) (80) 

A;^l fc— Tir) 
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where rip = 1 + [(J^ — si)/A],.aund> tit = 1 + [(m^ — si)/^]roundEIl ^ denotes the width of 
the mesh which is identified with the bin size in the data. The mesh points in the sums are 
determined by t^. =ti + {k- 1)A, k= 1,2..., with ti = 0.0125 and A = 0.025. The numerical 
coefficients gi^ take the values 3fc = 1 for 1 < fc < np and gi = gn^ = 0.5. The factors ru^ are 
determined by 

m = Wr{tk) if np < A: < n,-. 

r)(; = Q.5wr{tk) if k = Up or k = n-r. 

Using formula II80II and taking into account the definitions II67II and II68II . we calculate the 
required averages 

/ Tip rip— 1 np \ 

Gf = {Af Y. aK [fc] + 2 E E 9k9iCM , (81) 

\k = l fc^l l>k J 

I n-r n-r — l n-r \ 

G| = (Af E "tK W + 2 E E '^^'"Cfci , (82) 

\ fc — Tip k^rip Z — fc + 1 / 

GTG^ = (Af E E 9kViCki. (83) 

fc^l i — Tip 

We are now in a position to determine the uncertainties on the values of the "experimental" 
Adler function. They are induced from the errors of the experimental spectral function and 
from the errors on the parameters A and Sp. Let us represent again the "experimental" Adler 
function as a sum of the two terms showing explicitly the dependence of the terms on the 
parameters and on the spectral function 

£'"exp" (Q^ ^^ Sp :-yi) = Dcxp(Q^Sp :^'l) + DpQCD((^^vl^Sp), (84) 

the experimental and pQCD parts of the function arc determined as 

Dcxp (Q^ Sp : i>i ) = f " K.(Q^, t)vi {t)d t (85) 

Jo 

K.{Q^ ,t)vt^'^ (t)dt, (86) 

-p 

where IC{Q^,t) = 2Q^/(t + Q^f , vi{s) denotes the spectral function measured on the exper- 
iment (vi(s) = i'™''(s)) and «^^^(s) is the approximation to the spectral function evaluated 
within APT. Consider small deviations of the spectral function and the parameters from their 

mean values 

vi{t)=vi{t) + Svi{t), sp = sp + 5sp, A^=A^ + 5A^, (87) 

the change of the "experimental" Adler function under these variations is 

5D ..,,p. = 5Dexp + <5£>pQCD , (88) 

here we have used abbreviations Dcxp = Dcyip(Q^,Sp : vi) etc. The right hand side of II88II 
can be evaluated using formulas I I85II and I I86II . Preserving terms linear in the variations Svi, 
5sp and 5A^, we find 

<5i3"cxp" = <5„i Ocxp + E^^Ssp + Ej2SA^ (89) 

where 



Sv^D^^p = I " K{Q^ ,t)Svx(t)dt, 
Jo 



(90) 



Esp = K.{Q ,Sp)vi{sp)-\ — , (91) 

osp 

„ dPpQcujQ^A^Sp) 

^^' = dA^ • (^') 



^'^ here the subscript "round" refers to the integer nearest to the number inside the square 
bracket. 
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Using the trapezoidal rule, we approximate the integral on the right side of Eq. II90II 

rip 

S^,D,^p ^ G-i = ^ <?fc/C(Q2,tfe)M(ifc), (93) 

where the quantities rip, g^. and t^. are defined below formula II8OI I. To calculate the partial 
derivatives on the right hand sides of I I91I I and I I92I I. we use explicit formula II25I I for the pQCD 
part of the Adler function. We then obtain 



(Q2 + Sp)2^ 

1 r°^ 



agpQCD(Q^^^^p) _ 1 r e^PcSJt) ^^ 



where Pcff (*) = Pcff(^^£*) a-i^d to derive the last formula we have used the relation 

dr(sp) 1 

which can be easily derived from the definition II23I I. The mean squared deviation of the 
"experimental" Adler function is then determined as a sum of the six averages 



(5-D-exp" r = { S.,Dc.p r + E %{&sp)^ + E% (5vl2)2 

+2Esp£;^5sp5/l2 +2£;^p(5„il5cxp)5sp + 2£>p-(<5^iDcxp)<5/12. 

With the aid of formula II93I I. the first term on the right of Eq. II94II can easily be expressed in 
terms of the errors o-„j and the covariance matrix Cfe; 

( rip Tip — l Tip 

(<5„i£>oxp)2 = /i2j Y,glK\Q^,tk)cTl^[k] + 2 J2 Yl 9kgifC{Q\tk)JC{Q'',ti)Cki 

(_fc=l fe=l !=fe+l 

'(95) 

The second and third terms on the right of I I94II are determined in terms of the errors crgp and 
o'y^2 which we have already evaluated above (see II79II ). In order to evaluate last three terms 
on the right of 1194 l l. we use explicit expressions for the deviations 5sp and SA^ 

5sp=V-'(B2Gi-BiG2) 
M^ = ©- 1 (Ai G2 - A2 Gi ) , (96) 

the solution to the system II78II . This enable us to write 

<5sp<5yl2 = V-^iiBiAi + BiA2)GTG^ - B2A2GI - BiAiCf}, (97) 

the averages on the right hand side of II97II have been evaluated above (see Eqs. 1I8III . II82II 
and Il83l l). It remains to calculate the last two averages on the right hand side of I I94II . Using 
formulas II96I I we find 



(5„iDoxp)5sp = G3<5sp = ©-^(BaGiGs - B1G2G3) (98) 



(<5„iDexp)5yl2 = G3M2 = 0-i(AiG2G3 - A2G1G3), (99) 



employing now the trapezoidal sums I I8OI I and II93I I. we determine the averages G1G3 and G2G3 
in terms of the correlation coefficients Rj, ; 

Tip Tip 

GTG^ = A^J2 9k<^-i ['^l 12 9iK^{Q^,tink,i^v^ [I] 

k=l 1=1 

TlT "P 
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C Standard CIPT Consideration 

It is instructive to compare the modified procedure for extracting the coupling constant with 
the standard procedure formulated within conventional CIPT in the MS scheme. The t decay 
rate to the non-strange hadrons in the vector channel is given by |6] 

Rt,V= ^|Vud|'5Ew{l + 5QCD+<5Ew) (100) 

where (5qcd represents the QCD corrections, |Vmj| = 0.9746 ± 0.0006 is the flavor mixing 
matrix element, Sew = 1.0198 is an electro-weak correction term and (5ew ~ 0.001 is an 
additive electroweak correction (for these values see [9]). The QCD contribution is the sum 

5qcd=<5<°'+5(2)-h5np, (101) 

where (5''^' is the purely perturbative contribution, 5'^' is the dimension D = 2 effects from 
light quark masses, and (5np is the total non-perturbative contribution: (5np = <5('*' -t-i'^' -|-5(*' 
((5'^' are the OPE terms in powers of m^ ). We will use the estimates 5'^' = (— 3.3±3)xl0~'* 
and (5np = 0.0199 ± 0.0027, the ALEPH results obtained within the CIPT approach [9]. The 
experimental result for &^^' can be determined from the experimental spectral function via the 
relation 

1 + ^Si + 5(2) + 5np + 5ew = 4J°;'P, (102) 

where 



JTv — I Wt{s)v'{^'^ {s)d t 
Jo 



(103) 



and explicit expression of the function Wt{s) is given in II36I I. The relation II102I I follows from 
formulas II36II and lllOOII . Let us now determine the experimental error on J^^y induced from 
the experimental errors on i)™''(s). Using the trapezoidal rule, we replace the integral on the 
right side of Eq. II103II by the sum 

J^J ~AJ2 akwAsk)vi{sk) (104) 

fe=i 

where Af^ = 1 -|- [(m^ - si)/A],.ound, s^ = si -|- (A; - 1)A with si = 0.0125 and A = 0.025, and 
gj, are the numeric coefficients associated with the trapezoidal rule. From formula 11104 l l one 
easily evaluates the standard error on J'^^y 



-iKJ) 



y^ y^ gk9n.Wr{s^:)Wr{s„ 



(105) 



where C denotes the covariance matrix C^fe = E[{vi{si) — vi(si)){vi{sf^) — i'i(sfc)] which is 
available in [Toj. It follows from Eqs. IIIOIII and II102II that 

'^(Sqcd) = [a\5'') + a\5t,p)]'/''=4a{4J), (106) 

where we have ignored the small correlation between 5'") and Sj^p. With the data provided 
by ALEPH [10], from Eqs. |[T02t . itTOSll and HMl we obtain [3 



Si% = 0.2091 ± 0.0065cxp, (107) 

htly large value and error have been obtained, namely, S^Jp = 
0.2093 ± O.OOScxp. The perturbative QCD correction obtained within CIPT is represented via 



it should be noted that in |12| slightly large value and error have been obtained, namely, S^Jp 



^* Alternatively, we could have determined the error on <5cxp directly from the known error 
on Rr,v using formula UOOI I with the estimate Sj^p = 0.0199 ± 0.0027. 
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Table 14 Numerical values for the QCD scale parameter and strong coupling constant in the 
MS scheme for three active flavours extracted from the non-strange vector t lepton data within 
the conventional CIPT approach. The results obtained in consecutive orders of perturbation 
theory are given. The error bars refer to the experimental uncertainty only. 



Perturbative orders A GeV Qs{m^) 

LO 0.604 ±0.023 0.485 ± 0.019 

NLO 0.469 ±0.018 0.390 ± 0.011 

N^LO 0.430 ±0.016 0.367 ± 0.009 

N^LO 0.407 ±0.015 0.354 ± 0.008 

N*LO 0.395 ±0.015 0.347 ± 0.008 



the contour integral in the complex momentum squared plane |23II24| . This integral can be 
rewritten as 

4°' = ^ /" ^^ {(1 - ^'*') (1 + '='"')' -^RGlC^oe'^)} d<p, (108) 

where sq = m^ and djiQi{z) denote the RG improved perturbative correction to the Adler 
function defined in II14II . To calculate rfRGl(-z)i usually, the four-loop order RG equation is 
solved numerically for the running coupling. We find convenient to use the implicit solution to 
the RG equation at the four- loop order (relevant formulas can be found in [33! ) . The running 
coupling satisfies a transcendental equation which is solved numerically. To extract the value 
of the QCD scale parameter A = A-^^, one solves the equation 

5^^^{A) = si%. (109) 

Numerical values for the QCD scale parameter and strong coupling constant (for nf = 3 
number of flavours) extracted from the experimental value II107I I are given in Table 1141 We 
have used various approximations to the Adler function evaluated with the four-loop running 
coupling. For the unknown N^LO coefficient of the Adler function, we have used the geometric 
series estimate ds ~ 378 ±378 [12]. 
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